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We report numerical experiments of optical wave propagation in composites of high refractive 
index dielectric rods at frequencies where their first electric and magnetic Mie resonances are excited. 
The arrays of these particles have been extensively studied and proposed as non-absorbing and 
isotropic metamaterials. We show that negative refraction, observed in ordered particle arrays, is 
due to diffraction and that an effective medium theory yields constitutive parameters that do not 
reproduce the observations in these composites, whose transmission also depends on the sample 
shape. This is further confirmed by disordering the arrays, a case in which large transmission losses 
appear due to extinction by resonant scattering from the particles. Therefore, these composites 
although little absorbing have large extinction due to scattering. 
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INTRODUCTION 



Progress in metamaterial design, from microwaves to THz and optical frequencies \X~M^ based on the magnetic 
response of arrays of wires and split ring resonators (SRR) and their modifications to shorter electromagnetic waves, 
shows as essential limitations the existence of anisotropy and large absorption losses [g^. More recently, alternative 
structures based on Mie scattering by dielectric spheres [6-11] of relative large refractive index were extensively 
studied and proposed as models of lossless composites at microwaves; these were also extended to semiconductor 
cylinders and spheres, which were proven to possess similar resonant characteristics in the infrared and visible regions 
Although initially these latter structures relied on only the electric dipole or on multipole modes, later the 
possibility of exciting the first Mie magnetic dipole resonance was realized; which makes these particles equivalent 
from a fundamental point of view to the microwave wire and SRR elements as far as their electric and magnetic 
responses are concerned. Moreover, these cylinders and spheres provide isotropy 2D and 3D in addition to their 
resonance being subwavelength, both for the electric and magnetic excited dipoles. The purpose of this paper is to 
investigate the transmittance properties of media composed of these dielectric particles and in particular, whether 
they behave as uniform media at frequencies where their electric and magnetic dipoles are excited. Therefore this 
work is a test and assessment of whether the excitation of electric and magnetic resonances of high index Mie spheres 
and cyhnders constitutes an alternative with low losses and high refractive transmittivity, to previously developed 
metamaterial models. 

In this connection, we emphasize that in spite of the exhaustive studies on the possibilities of these resonant particles 
as metamaterial building blocks in and references therein, no such a test has been carried out. 

Although a unified view of the effective refractive index of metamaterials made of ordered arrays of such particles, 
as well as of photonic crystals (PC) producing negative refraction [14] was established fTE], the characterization of 
these composites as effective homogeneous media to the propagating wave has generally employed the method of 
the scattering (S) parameter by inversion of the complex transmittances and reflectances pL|]. Then, it remained the 
question of whether the effective constitutive parameters eg/ / and /ig/ / derived from effective medium homogenization 
procedures commonly employed [1, [sl, [U, [ill , were the same as those obtained by those methods that took into account 
the wave interaction with the microstructure of the composite unit cell. The answer, recently given [l6j, is negative. 
The linear dimension of such unit cells, which is typically between A/10 and X/6, is a much larger number than 
those of atom or molecule arrangements in transparent dielectrics acting to the pass of light. In fact, the maximum 
lattice constant versus wavelength for a metamaterial to behave as an uniform effective medium in negative refraction 
experiments, was established in [17|. Therefore, one may ask whether reducing absorption of the composite, like with 
these Mie resonances of dielectric spheres and cylinders, is sufficient to achieving an application as a metamaterial 
which may be considered as a refractive element. 

To this end, we shall see in this paper that, although contrary to metal elements like those of earlier left - handed 
material (LHM) designs, these dielectric cylinders and spheres grouped as "meta - atoms" of composites do not present 
absorption, their scattering cross section is quite large. Therefore, we will demonstrate that this time the losses of 
light transmission come from their large extinction cross section due to scattering, which moreover, when random 
arrangements of these particles are also addressed, yields a rather short transport mean free path, even if the medium 
is homogenized, and hence the transmittivity of a propagating beam through these media is low. 

As a 2-D equivalent to spheres, both ceramic [11] and of Silicon 0, rod array composites were suggested to exhibit 
metamaterial left-handed behaviour in the microwave and in the visible to mid-infrared (IR) ranges, respectively, due 
to the excitation of the cylinder magnetic Mie resonance. However, further studies pointed out [18] the necessity of 
more research to short out whether the backward wave behavior inside a periodic array of such Si cylinders is due to the 
band structure in the diffraction regime flA l or to a pure left-handed effect in the long wavelength range., even though 
experiments of microwaves in prisms of large permittivity ceramic rod arrays, either ordered or random, suggested a 
left-handed behavior [11]. 

The typical filling fraction / of the studied dielectric cylinder arrays |llMl3| is moderate / ~ 0.30, but the wavelength 
in vacuum to rod radius ratio: A/r was 61 in e. g. the experiment of [11] (A = 41.64mm), the lattice constant a to 
A ratio being 0.07, which is well below the aforementioned ratio a/ A = 0.1; whereas A/r = 9.8 and a/A = 0.45 in 
the model of the ordered Si cylinders of [lH (A = 1.55/im). On the other hand, the equivalent homogeneous media 
obtained from Snell law had an index of negative refraction n —0.6 for the PC slab of [12] and n —1.08 for the 
PC prism of [11], (A^ — 2.58/im for the case of [12], whereas A^ — 38.5mm for [H], A^ = A/n). None of these values 
is maintained when one changes the sample geometry, as we shall prove in this paper. 

In the homogenization procedure employed in the ceramic composite of [ll|, the effective parameters obtained in 
the band of left-handed behaviour have negative values both for the real part of the magnetic permeability, M^j, 
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as for that of the dielectric permittivity, e^j, with /i^j << e^j. However, near the resonance wavelength, the 
imaginary parts e^jj and /i^jj of both constitutive parameters are non - negligible compared to those real parts; 
this conveys high extinction by the composite material, which is not due to absorption losses as in previous models of 
LHMs composed of metallic elements, but produced by scattering from the high index particles of these metamaterials. 
The same happens in the near IR for the composite material of Si cylinders [12] and [13]. This will be further studied 
in this work. 

In this paper we carry out 2 - D numerical experiments with the finite element method (FEM) of propagation, of mid 
- IR waves, through a composite medium made of dielectric rods either in periodic or random positions. Nonetheless, 
for the sake of comprehensiveness, we shall also address the related problem of microwave propagation at larger scales. 
From the above discussion, our aim is to assess for the first time to what extent these structures may constitute a 
metamaterial model that overcomes the losses Him of previous composites. This conveys to discuss the validity of 
establishing effective constitutive parameters. The results should closely predict laboratory experimental observations 
because the calculations involved are exact. Since in this model the size of the samples are not huge compared to the 
wavelength [l6|, [l7], in order to assess whether the propagation depends or not on this size and on the shape of the 
sample, we employ two of such bodies for observing transmission: a rectangular block, or thick slab, and a prism, (i. 
e. in the 2-D calculations the latter being a triangle), composed of either a periodic or a disordered array of rods in 
air. Like in previous studies on these composites, the incident wave is assumed to be linearly polarized with electric 
vector E along the cylinder axes. 



SETS OF ORDERED AND RANDOM RODS AS METAMATERIALS 



Numerical procedure 

Maxwell equations are solved by using a finite element method (FE) (FEMLAB of COMSOL, http:/ /www. comsol.com). 
The calculation domain is meshed with element growth rate: 1.55, meshing curvature factor: 0.65, approximately. 
The geometrical resolution parameters consist of 25 points per boundary segment to take into account curved geome- 
tries in order to adapt the finite elements to the geometry and optimize the convergence of the solution. The final 
mesh contains about 10^ elements. To solve Helmholtz equation, the UMFPACK direct is employed. The boundary 
conditions of the simulation space are established both to keep the calculations from undesired window refiections and 
to avoid possible geometrical discontinuities. We then ensure that no inconsistencies due to properties discontinuities 
of the objets under study appear, and possible systematic errors are avoided. 

In the 2-D configuration, light depolarization is prevented by launching linearly polarized light with propagation 
vector in the XY - plane of the cylinder cross sections. Beam profiles are either, rectangular (plane waves): Eq exp(z(ki- 
T{ — cjt)), their widths being that of the simulation window, or Gaussian: Eq exp(— |R — Rop/2cr^) exp(z(ki • r — ujt)), 
r = (R, 2;). R = {x^y) and Rq = (xo,^o) are the transversal components of r and tq, respectively, a is the standard 
deviation or beam waist, and ki is the incident wave wavevector with |ki| = 27r/A. For the wavelength A, we shall 
address values either in the microwave or IR regions. The direction of propagation ki of such beams is thus normally 
incident to the OZ - axis of the infinite cylinders. The incident wave amplitude is normalized to |Eo| = IV/m (SI), 
which corresponds to a magnitude of the time average energy fiow | < S > | ~ 190W/m^. The criterium to choose 
the beam profile has been based on the major response of the composite characteristics to analyze: directionality 
of propagation in the media resulting from an homogenization method, and extinction processes when the inner 
structure of the rod distribution is considered. In the first case, we employ rectangular incident beams; otherwise, we 
use Gaussian beams. 

The results are thus expressed in terms of either the electric vector E(R) which points along the cylinder OZ - axis, 
the squared root of its time - averaged energy |E(R)|, the magnetic vector H(R) = {Hx^Hy) or the time average 
energy flow < S(R) >. These latter two vectors of course being both transversal, namely, in the XY - plane of the 
images to show next. 

Finally, to classify the whispering gallery modes (WGM) associated to the Mie resonances of the cylinders we will 
use the subscripts (i, j)^ i and j standing for their angular i - th and radial j - th orders, respectively. 



Electric and magnetic dipolar response of a dielectric cylinder in the microwave and mid - IR regimes 



We first address the response of one single ceramic cylinder of Bao^^Sro^^TiOs (BST) to microwaves. Figures 1(a) 
- (c) show the electric, magnetic and time - averaged Poynting vector distributions for linear polarized illumination 
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FIG. 1. (a) Electric field ^^(R) in a cylinder of BST ceramic with a dielectric permittivity e = 600 and radius r = 0.68mm. 
(b) Magnetic field H(R) (arrows) and its X - component (colors), (c) Averaged energy flow < S(R) > (arrows) and its norm 
(colors). In these figures, an s - polarized Gaussian beam of amplitude A — IV /m and standard deviation a — 12mm at 
A = 41.638mm is launched upwards (i. e. with Ri along the OY - axis), from below the cylinder, exciting its WGM : TMi,i. 



with the electric vector along the cylinder axis. Figure 1(a) exhibits Ezi^)^ indicating two electric currents flowing 
along the cylinder OZ - axis, one upwards and one downwards, respectively, and centered near opposite sides of the 
rod periphery, which corresponds to the E - spatial distribution of the dipolar WGM : TMi i. Figure 1(b) shows the 
magnetic vector lines in the XY - plane, characterized by arrows circulating around these electric currents, according 
to Ampere's law and behaves as that of a magnetic dipole. The time - averaged energy flow, characterized by the 
mean Poynting vector, is shown in Fig. 1(c) exhibiting an interesting circulation around the equatorial extremes of 
the cylinder section. 

On the other hand. Figs. 2(a) - 2(c) show the response of one single Si cylinder to infrared light. The characteristics 
quoted above for the BST ceramic cylinder are again reproduced for this Si rod, as expected from the analysis of their 
similar electric and magnetic resonances [l2j . 



Transmission characteristics of a slab with ordered rod distributions in the microwave regime 



We now consider a thick slab made of a periodic array of the BST rods in air whose response was studied in Figs. 
1(a) - (c), (see Figs. 3(a) and 3(b)). An s - polarized light beam is launched on this system at angle = 20°. The 
map of transmitted ^^^(R) field is shown in Fig. 3(a). A transmission angle ^ ~ 0° is observed within the sample. A 
resonant field distribution appears within the rods, being somewhat similar to that of Fig. 1(a), taking into account 
the frequency shifts due to the presence of neighbor rods. On the other hand. Fig. 3(b) shows both the time - 
averaged electric energy |£^2(R)| and energy flow < S(R) > in a detail of the block of rods. The latter showing the 



Bragg directions of propagation inside the crystal. Nevertheless, in the experiment of [11|] (cf. Figs. 3 of Ref. |ll| ) 
on a prism of an identical array, a refraction angle of about 20° was obtained instead, which would correspond to 
rieff — —1.08. A 75% of transmitted energy is lost in the rod slab of Figs. 3(a) and 3(b). 

However, homogenization from an effective medium theory employed like in [ll|, [l3|, EH, 22 1 leads us to estimate 
for this composite the effective homogeneous medium constitutive parameters = —36, e^jj = 14.7; /i^j = —0.9, 
/igjj =0.11 (see Fig. 4). Such a uniform medium produces an angle of refraction which is very small, similar to that 
of Fig. 3(a) but does not reproduce the inner structure of the wavefield. This is shown in Figs. 3(c) and 3(d) which 
display ^^(1^) and a detail of [^^(R)! and < S(R) > for such an effective medium. 
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FIG. 2. (a) Electric field ^^^^(R) in a cylinder of Si with a dielectric permittivity e = 12 and radius r = 158nm. (b) Magnetic 
field H(R) (arrows) and its X - component (colors), (c) Averaged energy flow < S(R) > (arrows) and its norm (colors). An s 
- polarized Gaussian beam of amplitude A = IV /m and standard deviation a — 2792nm at A = 1.55/im is launched upwards 
(i. e. with Ri along the OY direction), from below the cylinder, exciting its WGM : TMi,i. 



Now a 95% of losses appears in the transmitted energy emerging from the homogeneous medium of Figs. 3(c) and 
3(d). These results manifest the very different observed refraction depending on the sample on use: whether it is 
the composite or we deal with its EMT homogenization. In fact, although not shown here, we remark that we have 
reproduced the numerical results of [jjj for such an array if the sample is a prism which, as said above, yields a 
refraction completely different to that of Figs. 3(a) and 3(b). 

Hence, we infer the questionable matching of an EMT with the observations in such a regular array of rods. This 
is not surprising in the light of the results of [16] and [l7|, since the sample in these cases is not much larger than 
the inner structure, at difference with light transmission experiments in usual light refractive elements of nanoscopic 
molecular inner structure sizes. 



Microwave transmission in a slab with a random distribution of cylinders 



Different realizations, obtained by randomizing the array of Figs. 3(a) and 3(b) are now addressed. We next study 
a disordered distribution of BST rods as those employed in Section ??. Now, there is 85% of transmission losses 
through the slab of these disordered rods due to extinction produced by scattering. Figs. 5(a) and 5(b) show the field 
Ez(R) and a detail of both |£^(R)| and < S(R) >, respectively. A large amount of scattered light is lost both above 
and below the slab upper and lower low reflection boundaries. Again some few rods exhibit the excitation of Mie 
resonances on illumination as in Fig. 1(a). Figures 5(a) and 5(b) yield no clue of a forward transmission direction 
into the air at the exit of the sample. This result is obtained with both a beam and a plane wave. Although not 
shown here, we should state that we observed that averaging the field Ez(R) over many realizations of the random 
distribution of rods does not yield a forwardly transmitted beam, characterized by < Ezi^) > 13-221, distinguishable 
from the scattered light. 
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FIG. 3. (a) Electric field ^^(R) spatial distribution of a microwave propagating through a slab of an ordered array of BST 
rods as the one of Figs. 1(a) - (c) (e = 600, r = 0.68mm). The lattice constant is a = 3mm. An s - polarized Gaussian beam 
of amplitude A = IV /m and width a = 4a is launched on the PC from the left at wavelength A = 41.638mm = 7.2GHz) 
and at an incidence angle = 20° with the X - axis, (b) Electric field norm |£^2(R)| (colors) and < S(R) > (arrows) in a detail 
of the block upper left corner, (c) Same as in (a) in a slab occupied by a homogeneous medium whose effective parameters 
obtained from the EMT are Ce// = -36 + il4.7, fieff = -0.9 + iO.ll. (d) Electric field norm \Ez(R)\ (colors) and < S(R) > 
(arrows) in a detail of the block upper left corner. 




012345678 
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FIG. 4. Estimation of both the complex relative permittivity Ceff and permeability /Jeff for the BST rod array in Fig. 3(a) 
obtained the effective medium theory, e^ff (broken line with triangles. Blue online), e^ff (broken line with dots. Green online), 
/ifyy (full line. Red online) and /j^iff (broken line. Black online). The broken with dots vertical Une (violet online) at frequency 
ly — 7.2GHz indicates both Ceff and /Jeff values for the medium of Figs. 3(c) and 3(d): Ceff = —36, ^iff = 14.7, /i^ff = —0.9 
and /ilff = 0.11. The full vertical (blue online) line belongs to the jump of ef//. 
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FIG. 5. (a) Electric field ^^(R) spatial distribution in a random configuration of cylinders like those of Figs. 1(a) - (c) keeping 
the same filling fraction and the same conditions of illumination as in the ordered array of Figs. 3. (b) Electric field norm 
(R) I (colors) and < S(R) > (arrows) in a detail of the block central region which includes its upper transparent side. 



Ordered and random distributions of Silicon rods in the mid - infrared regime 

We next address a thick slab of ordered Si cylinders like that studied in Figs. 2. This arrangement is an extension 
to mid - IR ([12]) of the previous calculations at microwaves discussed in Sections ?? and ??. As shown in Figs. 
6(a) and 6(b) there appears negative refraction (i. e. a backward wave) in the block and no appreciable transmission 
losses. The positive and negative peak values of the wavefronts transmitted inside the array coincide with those 
cylinders that appear resonantly illuminated thus exhibiting the Mie Ti^i resonance like in Fig. 2(a). Mainly, there 
are < S(R) > arrows pointing in the direction normal to these wavefronts. Light is transmitted into the air side on 
the right through this slab with the same angle as that of incidence. This confirms the analysis of [12] for this crystal. 

On the other hand, Figs. 6(c) and 6(d) display the transmission of the same incident wave through a thick slab 
occupied by a uniform medium with a refractive index: n = — 0.36 + iO. 25/8 which resembles the propagation through 
the ordered array shown in Figs. 6(a) and 6(b). The real part = —0.36 has been estimated from Snell's law applied 
to Fig. 6(a), whereas the imaginary part = 0.25/8 was fitted to a transmittance of this homogeneous medium 
being approximated to that of Figs. 6(a) and 6(b). Notice that the transmittivity of the homogeneous slab of Figs. 
6(c) and 6(d) is smaller than that of the ordered array of Figs. 6(a) and 6(b). However cannot decrease much 



beyond 0.25/8, since otherwise well known instabilities due to divergences at the right side of the slab [28|-[30| appear. 
However, it should be remarked that this value of n does not coincide with that of an EMT, which according to the 
values of a/A and r/A for this array, as discussed in Section ??, are too large for a homogenization procedure to 
work with such a structure. This is further discussed next by employing other arrays of these cylinders with the same 
filling fraction and lattice parameter; namely, a thick slab of disordered Si rods and a prism of ordered, or disordered, 
cylinders. 

When these cylinders are disordered keeping / ^ 0.30, as seen in Figs. 7(a) - (c), there is no such negative refraction 
inside the block as that shown in Figs. 6. Now there is a huge extinction of energy due to high scattering by the 
particles. By averaging over several realizations, there is no observed refractive transmission of a forward component 
< Ez > into the air side. The block of disordered rods now scatters, see Figs. 7(a) and 7(b), more than 90% of 
the transmitted intensity both above and below its upper and lower low reflection boundaries. Hence, only some few 
spheres are illuminated showing their TMi^i Mie resonances. Figure 7(c) shows that | < S(R) > | is transmitted into 
the air region at the right of the slab with less than a 1/10 of its incident wave value. 



A prism of ordered or disordered Si distributions in the infrared 

Let us consider now a prism of a uniform medium of refractive index n = —0.36 + iO.25/8. Notice that this value of 
n is the one of Figs. 6(c) and 6(d), derived from refraction at the ordered array of Figs. 6(a) and 6(b). As shown 
in Figs. 8, refraction into the air takes place with Oi = 18.4° and Of = 38.07°, the latter being the negative angle of 
refraction at the prism larger side. This is shown in Figs. 8(a) - (c) and serves us as a reference to study transmission 
through both an ordered and disordered array of Si rods contained in a sample with this prism geometry, which also 
was the one employed in transmission observations at microwaves in [Tlj . 
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FIG. 6. (a) Electric field ^^(1^) propagating in a thick slab occupied by an ordered array of Si rods (e = 12 and r = 158nm; the 
lattice constant is a = 698nm). An s - polarized Gaussian beam of amplitude A = IV /m, a = 4a and wavelength A = 1.55/im 
is launched on the slab from the left at an incidence angle = 20°. (b) Electric field norm |^;z(R)| (colors) and < S(R) > 
(arrows) in a detail of the block upper left corner, (c) Electric field ^^(R) in a uniform medium optically equivalent to that 
of (a) and (b), whose electric permittivity is Ceff = —0.36 + iO.25/8. The real part of such a value is estimated by Snell's law, 
(Oi — 21.54°, Ot — 38.07°). The imaginary part has been estimated from the transmittivity of Figs. 6(a) and 6(b). (d) Electric 
field norm |£^2(R)| (colors) and < S(R) > (arrows) in a detail of the block upper left corner. 



Figures 9(a) - (c) show the wave propagation on illumination of the prism filled with an ordered array of Si cylinders 
identical to that of the thick slab of Figs. 6(a) and 6(b). As shown, there is now absence of negative refraction at the 
larger side interface of this prism. This contrasts with the observation in the sample with the same shape filled with 
a uniform medium, as displayed in Figs. 8(a) - (c) and also with the case of a thick slab filled with the same array, 
as seen in Figs. 6(a) and 6(b). Hence, at difference with the crystal of Si rods in the block, the same crystal in the 
prism does not reproduce negative refraction, but rather a set of diffracted orders into the air which are associated 
to the prism angle a, according to the conservation of the transversal wavevectors at the larger interface. Also, the 
wave propagation inside the prism, which should be like that appearing in the thick slab of Fig. 6(a), but now at 
the same direction as the incident wave, (since now this latter wave incides on the prism at Oi = 0°), is observed in 
Figs. 9(a) and 9(b) to be quite different, with a complicated structure due to the interference of the different Bragg 
waves. This result, once again, points out the dependence of wave propagation not only on the inner structure of the 
composite, but also on its shape. In addition, these observations demonstrate that the negative refraction found in 
the slab of ordered Si rods, (cf. Fig. 6(a)) is not due to an effective homogeneous medium effect, but it rather comes 
from a consequence of the diffraction in the array of Si rods due to its lattice symmetry. This answers the question 
posed in [18]. 

This latter remark is further confirmed by disordering this Si cylinder array within the prism, (see Figs. 10(a) - 
(c)). Again, no refracted beam into the air is observed, and a large portion of energy is lost by scattering from the 
composite elements. Many outgoing beams mix with each other, so that no forwardly transmitted beam inside the 
prism, characterized by < ^^^(R) > can be distinguished by averaging over several realizations of the random array. 
The apparent refracted beam at Ot = 0° shown in the air side exiting the prism, changes as one varies the random 
realization of rods. 



9 



(a) 




\p o Oq ° O <^ O c 




o '5 <^ o o o v-o 

n a£ <= „ ° O C O Son 
nO n° ° C C, 




o c o 0^° o « c o ° o^o 
j>cr' o o o o" c 







FIG. 7. (a) Electric field ^^(R) when the configuration of Si rods shown in Figs. 6(a) and 6(b) is disordered, keeping the same 
filling fraction / = 0.30 and the same conditions of illumination, (b) Electric field norm |^2(R)| (colors) and < S(R) > (arrows) 
in an inset of the block left region which includes its left interface, (c) Map of the averaged energy flow norm | < S(R) > |. 



CONCLUSIONS 



In this paper we have assessed the transmittance of composites of dielectric particles whose Mie resonance electric 
and magnetic modes were proposed by previous extensive studies to yield negative refraction. We have shown that 
such structures cannot be homogenized neither reproduce the propagation observed in the frequency regions of study. 
In this way, we have proved that the negative refraction previously found in ordered arrays is thus a diffraction effect 
which disappears as soon as the particle distribution is randomized, and does not reproduce the transmission of their 
corresponding EMT uniform media, this was an open question so far. 

In addition, we conclude, first, that the effective parameters obtained from a homogenization theory, do not repro- 
duce the propagation observed through an ordered array of these elements. This is further seen when these arrays 
are randomized. If an EMT worked, it should not depend on whether the "meta - atom" distributions were ordered 
or disordered. Then strong scattering by disordered rods extinguishes most of the incident energy, and there is not 
negatively refracted forward beam observed. 

Second, the behavior of the transmitted beam in ordered arrays also depends on the shape of the sample and hence 
does not match with an EMT which does not include this shape. These conclusions are consistent with the well 
known difficulty of working out EMTs within the frequency range of resonance of the composite elements [igI] . Similar 
consequences should hold for 3D composites of resonant dielectric spheres. 
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FIG. 8. (a) Electric field ^^(R) in a prism (top angle a = 18.4°) occupied by a uniform medium whose electric permittivity 
is Ceff = —0.36 + iO.25/8 illuminated from the left. This prism is optically as dense as the homogeneous thick slab of Figs. 
6(c) and 6(d). (b) Electric field norm |^2(R)| (colors) and < S(R) > (arrows) in a detail of the prism upper region, (c) Map 
of the averaged energy flow norm | < S(R) > |. An s - polarized Gaussian beam of amplitude A = IV /m, a = 4 x 698nm and 
wavelength A = 1.55/im is launched from the left on the prism at Oi — 0° with its left side. 
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FIG. 10. (a) Electric field ^^^(R) in the configuration shown in Figs. 9(a) - (c), (now with a top angle a = 22°), the Si 
rod array being now randomized but keeping the same filling fraction / = 0.30 as in the ordered array. The conditions of 
illumination are the same as in Figs. 8(a) - 8(c). (b) Electric field norm |^2(R)| (colors) and < S(R) > (arrows) in a detail of 
the prism upper area, (c) Averaged energy flow norm | < S(R) > |. 
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